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Abstract. - We discuss the superfluid phase transition of a strongly interacting Fermi gas with 
unequal (asymmetric) chemical potentials in two pairing hyperfine states, and map out its phase 
diagram near the BCS-BEC crossover. Our approach includes the fluctuation contributions of 
"preformed Cooper pairs" to the thermodynamic potential at finite temperature. We show 
that, below a critical difference in chemical potentials between species, a normal gas is unstable 
towards the formation of either a finite-momentum paired Fulde-Ferrell-Larkin-Ovchinnikov 
superconducting phase or a uniform superfluid, depending on the asymmetry and interaction 
strengths. We determine the value of critical chemical potential mismatch, and find that it is 
consistent with a recent measurement by Zwierlein et al. [Science 311, 492 (2006)]. 



Recent experimental advances in manipulating ultracold atomic Fermi gases via a Feshbach 
resonance have attracted a great amount of interest in widely varying fields from condensed 
matter physics, atomic molecular and optical physics, to particle and astro physics [1-3]. 
Thanks to this precisely controllable environment, it is now possible to experimentally ex- 
plore fundamental problems of many-body physics in the strongly interacting region. An 
example of particular interest is the understanding of fermionic pairings of two-component 
Fermi gases with mismatched Fermi surfaces or chemical potentials, in the crossover from 
the weak coupling Bardeen-Cooper-Schrieffer (BCS) regime to the strongly correlated Bose- 
Einstein condensate (BEC) regime [4-12]. 

In standard BCS superconductivity, the presence of a Fermi surface mismatch suppresses 
superfluidity, and the pairing mechanism may be qualitatively altered. Various competing 
states have been proposed within mean-field theory to describe the ground state of an asym- 
metric Fermi gas in the weak coupling BCS regime: the Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state [4], where Cooper pairs possess a finite center-of-mass momentum, the breached 
pair phase [5] and the Sarma phase [6,7], or phase separation phases with a paired BCS super- 
fluid being surrounded by an unpaired normal gas [8] . When the interaction strength is tuned 
across the BCS-BEC crossover, a rich phase diagram consisting of several these scenarios has 
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also been put forward [9-11]. Very recently, there have been two experimental investigations 
of fermionic superfluidity of ^Li atoms with mismatched Fermi surfaces [2,3]. 

In contrast to these earlier studies, in this paper we tackle the crossover problem of a ho- 
mogeneous asymmetric Fermi gases at finite temperature, starting from a well-defined normal 
state. Our goal is to determine the superfiuid transition temperature and to obtain a reliable 
phase diagram for these systems beyond mean-field. Our calculations reveal several notable 
features: (i) The normal state is stable above a critical chemical potential mismatch. The 
critical value predicted is in qualitative agreement with the recent measurement, (ii) On the 
BCS side, the critical value is roughly proportional to the transition temperature of a sym- 
metric gas and, therefore is exponentially small. Below it, the normal gas is unstable towards 
the formation of FFLO states. As the mismatch of Fermi surfaces decreases, a uniform BCS 
superfiuid is more favorable, (iii) In the strong coupling BEG regime, the critical difference 
in chemical potentials is of the order of binding energy. The superfiuid state, well described 
by a mixture of tightly bounded Cooper pairs and unpaired fermions in this limit, is thereby 
remarkable robust, (iv) The transition temperatures at fixed chemical potential imbalances 
are also determined, giving rise to a finite-temperature phase diagram. 

All these results are derived below by generalizing a thermodynamic approach by Nozieres 
and Schmitt-Rink (NSR) [13-15] to asymmetric Fermi gases. This approach takes into account 
the large fiuctuation effects that are necessary in order to capture the essential physics at the 
crossover. 

NSR approach to an asymmetric Fermi gas. — Because the Feshbach resonance of ^Li atoms 
used in the experiment is extremely broad [16, 17], we can use a single-channel model to 
describe the Fermi gas across a Feshbach resonance: 

W = ^ (ek - Ha) cJ<,Cka + UY^ C+ C+_k^Cq_k'iCk'T- (1) 

kff kfc'q 

Here c^^ is the creation operator for the fermionic atoms, and the pseudospins <t =t, J, denote 

the two hypcrfinc states of ^Li. The masses are the same, so ek = ?i^k^/2m for both species, 
but the chemical potentials are different, i.e., fXf^i = M i "^A*) to account for the asymmetry. 
We will focus on the situation with fixed total particle number n = n-i + ni and fixed chem- 
ical potential difference. U is the effective interaction strength and is related to the s-wave 
scattering length a via the regularization: (Anh^a/m)^^ = U~^+ X^kl^^k)"^- A two-body 
bound state arises in vacuum once the scattering length a becomes positive. 

Within the NSR approach [13], there are two essential ingredients in the thermodynamic 
potential at a temperature T: Q{T,iJ.,dn) = + ^pf- Here fto corresponds to a free Fermi 
gas, 

no = ^J2^nf[- (ek - Mt)] + In / [- (ek - Mi)] , (2) 

k 

where /3 = l/fc^T and / (x) = [exp(/3a;) + 1]~^ is the Fermi distribution function, while ilpf 
is associated with the pairing fluctuation contributions, and may be determined by summing 
an infinite series of ladder diagrams [13] , 

-|-oc 

= -E / ^9{u;)d{ci,u;). (4) 

^ — OO 



XiA-Jl LlU^ and Hui Hu^'^: BCS-BEC CROSSOVER IN AN ASYMMETRIC TWO-COMPONENT FeRMI GAS3 



Here q is the center-of-mass momentum, — 27rn//3 is the bosonic Matsubara frequency, 
g (to) = [exp(/3w) — 1]^^ is the Bose distribution function, and the particle-particle vertex 



function r -'^(q, = 1/C/ 



may be written as 



E 



/„ - 2eq/2 - 2ek + 2^ 2ek 



(5) 



with — eq/2±k ^ T ^M- Following NSR, in Eq. we have converted the sum over if„ 
into a contour integral, and have rewritten flpf in terms of a phase shift defined by d{q,Lu) = 
— Iniln[— r^^(q, w„ to + iO^)]. The chemical potential n for our model is determined from 
the number identity n — —dfl/d^, or, 



n = n% (T, fi, Sfi) + 2nB (T, ^, S^i) 



(6) 



where uPp = T^i.afi'^^ - fJ-^), and ub = -d^pj / d{2fj.) = I/tt X^q /.^^ rf^g(w)9(5(q,w)/(9(2^) 
may be interpreted as the number of "preformed Cooper pairs" . The superfluid phase transi- 
tion occurs when the particle-particle vertex function develops a pole at w„ = for a certain 
value of q. Therefore, the transition temperature Tc can be conveniently obtained by the 
Thouless criterion, 

maxF^^ (q, Wn — 0) |t=Tc= 0, 



(7) 



which is generalized here to take into account a nonzero center-of-mass momentum. We must 
solve Eq. (0 together with Eq.Q self-consistently, to obtain fJ,{Tc) and Tc for given chemical 
potential imbalance S/j, and interaction coupling, where kp = (37r^n)^/^ is the Fermi wave 
vector. 

BCS and BEC limits. — The NSR formalism presented above is a simplified description 
of the full crossover problem for an asymmetric Fermi gas. It is asymptotically exact in the 
extreme weak or strong coupling limit {l/kpo, — > ±oo). In between, it is believed to provide 
a qualitative interpolation scheme [13]. 

In the weak coupling limit, the phase shift (5(q,CL;) is small: the number equation © 
reduces to n = rPp (T, /i, 5^). The chemical potential is therefore that of a non-interacting gas 
of fermions, /x ~ = h'^kp/2m. Keeping in mind that Sfx, Tc <C ji, to leading order of 1/(/3/j,) 
and qvp/^n, one finds the following expression for the inverse of two-particle vertex function: 



2^ 

mkp 



T 

\n— + hi (Sfi) 



2„2 



7i2 {S^i) , 



(8) 



where kgTc = k^Tgcs — {%/iT)e^^'^tpcyi\i{iT /2kpa) is the BCS transition temperature for 
a respective symmetric gas, hi {x) = — dy ln{l3y)[f' {y -I- x) -I- f'{y — x)] + j + ln(4/7r) 
and /i2 (x) = — dy{l/y)[f"{y -I- x) -I- f"{y — x)]. At small chemical potential imbalance 
{5(i/Tc 0), /i2 i^lj) > 0, and the superfluid instability occurs at q = 0. Asymptotically, we 
obtain 



Tc(M 
TP 



1 



7 
12 



8C' (-2) 



C(3) 



\kBT0j 



(9) 



where is the Riemann zeta function. As dfi increases, /12 crosses zero and then 
becomes negative. Hence the maximum value of F"^ locates at a nonzero center-of-mass 
momentum, triggering the emergence of a non-uniform FFLO state. The condition /12 i^l^) — 
defines a tricritical point connecting both two superfluid phases and the normal phase. We 
find numerically that S^tri — l-OTkgTc and Tc^tri — 0.56T^, in agreement with a previous 
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study from Ginzburg-Landau theory [18]. As Sfj, increases further and approaches a critical 
imbalance, the transition temperature shrinks to zero, marking a quantum phase transition 
to the normal state. In this case the expression ^ becomes inapplicable, since the resulting 
qvp may become comparable with In the limit of /3 — > oo, alternatively we may cast Eq. 
Ipl in the form. 



2*2 



rnkp 



m — h F ix) 



TT 



kBT 
5^1 



(10) 



where = 1 - 7 - ln(2/7r) - l/21n|a;2 - l| - l/(2a;)ln|(a; + l)/{x - 1)| and x = qvp/SfJ- 

The function F (x) has a maximum at x — Xc — 1.20, with F {xc) — 0.286. Thus the 
quantum phase transition takes place at S^c — ISSlksT^ ■ Slightly below S^c the transition 
temperature shows a square root behavior. 



1/2 



(11) 



Deep within the BEG regime, the results are greatly affected by the phase shift contribu- 
tions. The particle-particle vertex function develops a discrete pole at a positive frequency 
Lu = eq/2 — 2/1 + €{,, where = It? /mo? is the binding energy. Thus, all fermions in the less 
populated hyperfine state will pair up with atoms in other state to form bound pairs. As 
a result, in the low energy regime the phase shift is dominated by a bound state part, i.e., 
S{q,ijj) ~ TT0{eq/2 — 2/i + £{,), and the number of pairs consequently acquires a standard form 
of an ideal Bose gas: ub — ^ J^q ^(^q/^ ~2fi + £{,). Accordingly the chemical potential ^ 
asymptotes to half the binding energy. On the other hand, to accommodate the remaining free 
excess fermions with density Sn — n^i — > 0, the chemical potential imbalance S/i should 
be necessary large. One may anticipate that ^ — while /i| ~ epiSn), where epiSn) is 
the respective Fermi level of these excess fermions: 6n w /(^k ~ M ~ (^/i). At the transition 
point the Thouless criterion gives rise to fJ-(Tc) = —eh/2. The superfluid transition tempera- 
ture is therefore determined by the strong coupling number equation describing a mixture of 
non-interacting bosons and fermions. 



e('k-A'T)/fcBr<= -I- 1 ^2 



mkBTc 



3/2 



(12) 



By setting Tc = 0, we obtain a critical chemical potential imbalance: S^c = £6/2 -1- 2^^^eF. 

Two striking features emerge from this analysis: first, in both BCS and BEC-like phases, 
the superfluid state is destroyed by a sufficient large chemical potential imbalance. On the 
BGS side, the critical difference in chemical potentials is exponentially small, while deep 
within the BEG regime, it is set by the binding energy. Second, in the BGS-likc phase, the 
preformed Cooper pairs can adjust its center-of-mass momentum in response to the Fermi 
surface mismatch. This results an inhomogeneous FFLO superfluid state in the vicinity of 
critical chemical potential difference. With these background, let us turn to the numerical 
NSR calculations. 

Transition temperatures and phase diagrams. — Figure 1 presents our results for the su- 
perfluid transition temperature as a function of the chemical potential imbalance at various 
interaction strengths. The dashed lines shows the region with the onset of FFLO phases. 
We plot the number of "preformed Gooper pairs" in dotted lines. Two mechanisms for the 
depression of transition temperatures with increasing imbalance may be identified. On the 
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Fig. 1 - (color online) Superfluid transition temperature Tc as a function of the chemical potential 
imbalance at various interaction strengths: (a) l/kpa — —1, (b) the unitary limit l/kpa = 0, and (c) 
l/kpa = +1. The phase transitions into an inhomogeneous FFLO state have been marked by dashed 
lines. The dot-dashed lines in (a) shows, respectively, the asymptotic behavior near Sfi — and 
Sfi = Sfj,c in Eqs. Q and lllL while that in (c) shows the Tc determined from the strong coupling 
number equation H12|l . The dotted lines display the number of "preformed Cooper pairs". In (c) 
Tbec = (?i^/m)7r[n/(2<;(3/2))]^''^/fc_B is the transition temperature for an ideal Bose gas. 



strong coupling regime, the decrease of is accompanied with the reduction of number of 
tightly bounded pairs, as expected from the general picture of BEC. In contrast, in the oppo- 
site BCS-like phase, the number of pairs keeps almost constant. These pairs, resulting from 
the many-body effects in this case, would be very fragile with respect to the Pauli blocking 
effects arising from mismatched Fermi surfaces. We therefore explain the reduction of Tc as 
due to the loss of phase coherence between pairs. 

Figure 2a gives the critical chemical potential imbalance Sjic throughout the BCS-BEC 
crossover, constituting a zero-temperature phase diagram. For comparison, we plot the ex- 
perimental data on the critical Fermi energy difference 5Ep/ep. These data are obtained 
indirectly by assuming a non-interacting dispersion: 5Ep/ep = [(1 + 5c)^^''^ — (1 — 5c)^^^]/2, 
where 5c is the measured critical population imbalance (see, i.e., the Fig. 5 in Ref. [2]). We 
find that our predictions agree qualitatively with the experimental results in the BCS-unitarity 
regime. However, on the BEC side they are not consistent. To understand this discrepancy, 
several remarks should be in order: First, our calculation is for a homogeneous gas, while the 
experiment is done in a trap. The presence of the trap tends to yield a phase separation, 
which may further complicate the comparison. Secondly, the experiment is performed at a 
unknown temperature. A close examination of the experimental data in the unitary limit 
suggests that there is an appreciable effect due to the finite temperature. Finally, the ex- 
perimental data refer to the critical Fermi energy difference, instead of the critical chemical 
potential difference as we calculated. Only in the weakly coupling BCS regime do the chemical 
potentials equal the Fermi energies. Thus it appears to be the most serious reason responsible 
for the discrepancy between our predictions and experimental data on the BEC side. In this 
regard, a theoretical calculation for the critical population imbalance will be useful, enabling 
a direct comparison with the experiment. Unfortunately, in a region around the unitary limit 
(^0.5 < l/kpa < +0.2, see Fig. 2b), we find that the NSR approach generally leads to a nega- 
tive population imbalance at a positive chemical potential difference, implying an un-physical 
compressibility, i.e., dSn/dSfj, < 0. This suggests the breakdown of the NSR treatment around 
the unitarity regime. We note that the negative compressibility is in close connection with 
the non-monotonous behaviour of the transition temperature predicted by NSR scheme for 
a symmetric Fermi gas [f3, 14]. Presumably, the breakdown of NSR approach is due to the 
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Fig. 2 - (color online) (a) Predicted critical chemical potential imbalance (solid line) in the BCS- 
BEC crossover. Solid symbols are the experimental data on the critical Fermi energy difference. 
The open circle with an error bar in the unitary is a rough estimate from quantum Monte Carlo 
simulations in Ref. [12]. The position of the tricritical points {5^tri) has been plotted by a dashed line. 
Inset shows the critical imbalance in the logarithmic scale, where two dotted-dashed lines represent, 
respectively, the critical imbalance in the extreme BCS and BEC limits, i.e., Sfic — LSSlksTBCS and 
6fic ~ et/2 + (b) Predicted critical population imbalance Sc . The NSR approach (at T = 0) 

breaks down around the unitarity regime —0.5 < l/kpa < +0.2, and thereby 5c is unknown. The 
dotted-dashed line shows the perturbation result in the BCS limit; Sc = 3/[2(l — 2kFa/TT)]Sfj.c. 



exclusion of the interactions between Cooper pairs, and may be avoided by a self-consistent 
improvement, with which the un-physical peak structure of Tc is shown to be removed [17,19]. 
In Fig. 2b, we report the critical population imbalance outside of the imitarity regime. The 
two experimental points on the BEC side now become consistent with our theory. 

In Fig. 2a we have also shown the chemical potential at the tricritical point (with nonzero 
Tc) by a dashed line. Therefore, the area enclosed by the solid line and dashed line marks 
the possible region for forming the FFLO state. The region ends up around the unitary limit. 
This fact, together with the low Tc for FFLO phases, may be used to understand the reason 
of why such states are not observed in the experiment. 

For completeness, in Fig. 3 we consider the transition temperature against the interaction 
strength at several fixed values of chemical potential imbalance, which makes up a finite 
temperature phase diagram. In particular, there is no superfiuid state on the BCS side above 

~ O.Sej?. 

Conclusions. — In summary, based on a generalized NSR approximation beyond mean- 
field, we have qualitatively determined the superfiuid transition temperature and a phase 
diagram for a two-component Fermi gas near a Feshbach resonance with mismatched Fermi 
surfaces. The resulting critical imbalance in chemical potentials agrees qualitatively well with 
the recent experimental findings [2]. While our approach starts from a well-defined normal 
state and thereby avoids the complication of assuming a possible pairing scheme, it cannot 
be used to identify the transition between different competing phases within the superfiuid 
domain. This fascinating issue could be addressed by extending the current NSR approach to 
the broken-symmetry state [20] . 

* * * 
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Fig. 3 - (color online) Finite-temperature phase diagram with the superfluid transition temperature 
at the crossover plotted versus chemical potential imbalances. For a better illustration, we have 
subtracted half the binding energy from the imbalance, except for the curve with 5fi = . 
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